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We study front propagation when an invading species competes with a resident; we assume nearest- 
neighbor preemptive competition for resources in an individual-based, two-dimensional lattice model. 
The asymptotic front velocity exhibits power-law dependence on the difference between the two 
species' clonal propagation rates (key ecological parameters) . The mean- field approximation behaves 
similarly, but the power law's exponent slightly differs from the individual-based model's result. 
We also study roughening of the front, using the framework of non- equilibrium interface growth. 
Our analysis indicates that initially flat, linear invading fronts exhibit Kardar-Parisi-Zhang (KPZ) 
roughening in one transverse dimension. Further, this finding implies, and is also confirmed by 
simulations, that the temporal correction to the asymptotic front velocity is of 0(t~ 2 ^ 3 ). 
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I. INTRODUCTION 

The dynamics of propagating fronts integrates con- 
cepts shaping our understanding of how invasive species, 
emerging infectious disease, and evolutionary adapta- 
tions spread across ecological landscapes 1]. Indeed, ob- 
jects as seemingly dissimilar as chemical reaction fronts 
0| and the spread of opinions share certain basic 
spatio-temporal properties. Fisher £| and Kolmogorov 
et al. first addressed velocity characteristics of simple 
fronts by modeling a favored mutation's one-dimensional 
spread with a reaction-diffusion equation. A lengthy se- 
ries of biologically generalized reaction-diffusion models 
has since appeared 0. However, recent developments 
emphasize the ecological realism of discrete individuals 
000 Our study analyzes front p ropaga tion when two 
plant species compete preemptively [Hjlllllllillilllil f° r 
a common limiting resource. Discrete individuals of each 
species propagate clonally, so that competitive interac- 
tions are spatially localized. An "invader" species has an 
individual-level reproductive advantage over a "resident" 
species, so that competition is asymmetric. 

This paper focuses on one-dimensional fronts separat- 
ing the species in a two-dimensional environment. We 
study the asymptotic front velocity, as well as the tem- 
poral and finite-size corrections (or rates of convergence) 
to this velocity. Furthermore, we investigate roughen- 
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ing of the moving fronts from a non-equilibrium interface 
viewpoint 0,0]- Asymptotic properties of initially flat, 
linear fronts do offer insights concerning the competitive 
dynamics of locally propagating plants. Consider trem- 
bling aspen (Populus tremuloides) . Significant seed pro- 
duction, hence long-distance dispersal, occurs only about 
once each five years Most growth is clonal, where 

a new tree grows from an existing tree's roots. Clonally 
propagated trembling aspen clusters occasionally expand 
to several thousand trees, and cover > 40 ha To 
model such systems, one can assume that introductions of 
an invader by seed occur rarely and stochastically, in both 
space and time. We have shown 0, EJ E, that the 
time evolution of the invader and resident populations in 
such systems can be well described within the framework 
of homogeneous nucleation and growth |2l| . In particu- 
lar, in two dimensions, for sufficiently large systems, the 
typical time (lifetime) until competitive exclusion of the 
weaker competitor scales as t^/i; 2 ) -1 / 3 [§. IsL Il9l l2fj| . 
where / is the stochastic nucleation rate per unit area of 
the successful clusters of the better competitor, and v is 
the asymptotic radial velocity of the growing (on aver- 
age) circularly symmetric fronts. It is, thus, clear that 
the full understanding of the dependence of the lifetime 
on the local rates of the systems requires the knowledge 
of the velocity of the front separating the two species. 
Furthermore, as circular fronts grow sufficiently large, so 
that curvature corrections become negligible, frontal ve- 
locities approach values for linear fronts . 

Recently we have reported preliminary results on the 
front velocity in the model studied here |2^|. This paper 
extends our analysis by investigating the front's propaga- 
tion as a non-equilibrium interface |l5l Utij. Our Monte 
Carlo simulations not only provide numerical estimates 
for the front velocity, but, through a detailed finite-size 
analysis, also identify the universality class of the mov- 
ing and roughening interface 0, 0] separating the two 



species. Our results indicate that the asymptotic front 
velocity scales as a power law with difference between the 
two species' clonal propagation rates. Further, we find 
that initially flat, linear invading fronts exhibit Kardar- 
Parisi-Zhang (KPZ) |'23| roughening in one transverse di- 
mension. This finding implies, and was also confirmed 
by our simulations, that the temporal correction to the 
asymptotic front velocity is of 0(t~ 2 ^ 3 ). 

We organize the remainder of the paper as follows. In 
Sec. II we define the spatially explicit, individual-based 
model of two-species competition. In Sec. Ill we compare 
simulation results for the asymptotic front velocity with 
results from the mean- field equations. In Sec. IV we carry 
out the analysis of the interface roughening characteris- 
tics of the front, also yielding the temporal and finite-size 
corrections to the asymptotic front velocity. We discuss 
and summarize our results in Sec. V. 



II. TWO-SPECIES INVASION MODEL WITH 
PREEMPTIVE COMPETITION 

Our analysis treats the velocity and roughening of in- 
vading fronts as functions of the propagation and mor- 
tality rates of each species, with possible temporal and 
finite-size corrections. On a two-dimensional L x xL y lat- 
tice, a site represents the minimal level of resources nec- 
essary to maintain a single individual. Competition for 
the resource is preemptive that is, a currently occu- 
pied site cannot be colonized by any species until mortal- 
ity of the occupant opens the site. The local occupation 
number at site x is n;(x) = 0,1 with i = 1,2, repre- 
senting the number of resident and invader individuals, 
respectively. Since two species cannot simultaneously oc- 
cupy the same site, the excluded volume constraint yields 
ni(x)ri2(x) = 0. A species may occupy new sites only 
through local clonal propagation. Therefore, a species 
occupying site x may only reproduce if one or more of its 
neighboring sites is empty (here we consider only nearest- 
neighbor interactions). 

During our time unit, one Monte Carlo step per site 
[MCSS], L x L y sites are chosen at random for updating. 
The local configuration of a chosen site is updated ac- 
cording to the following transition rates. An empty site 
may be occupied by species i through clonal reproduc- 
tion from a neighboring site at rate a^^x), with a,i be- 
ing the individual-level reproduction rate for species i, 
and r?i(x) = (I/4)E x / enn ( x )n i (x / ) is the density of species 
i in the neighborhood around site x; nn(x) is the set 
of nearest neighbors of site x. An occupied site opens 
through mortality of the individual; the mortality rate p, 
is the same for both species. The transition rules for an 
arbitrary site x can be summarized as follows: 

ai -H x) l, Q2 H X) 2, 1-^0, 2^0, (1) 

where 0, I, 2 indicates whether a site is open, resident- 
occupied, or invader occupied, respectively. 



One should note that the above discrete stochastic 
model, governed by preemptive competition, is a two- 
species generalization of the basic contact process 
13 13- Each species, in the absence of the other 
one, becomes extinct (through a transition to an absorb- 
ing phase |2f|) if cti<a c (p), where a c (p) « 1.65/Lt 
[and a c (p) = p in the mean-field approximation, see Sec. 
III. A]. We are interested in the scenario where compe- 
tition between the two species drives the dynamics (i.e., 
not extinction by insufficient colonization rates), and one 
of the species (the invader) has a reproductive advantage 
over the resident. Hence, we study the a c {p) < ot\ < eti 
regime. 

To study front propagation, we impose periodic bound- 
ary conditions along the y-direction of an L x xL y lattice. 
The initial condition is a flat linear front (straight ver- 
tical line), i.e., the invader completely occupies a few 
vertical columns at the left edge of the lattice, and all 
other remaining sites are occupied by the resident species. 
The direction of propagation, therefore, is along the x- 
direction. As the simulation begins, many individuals 
of both species die in a few time steps, and the den- 
sities on both sides of the front quickly relax to their 
"quasi-equilibrium" values, where clonal propagation is 
balanced by mortality [Fig.^. As the simulation evolves, 
we track the location of the invading front by defining the 
edge as the location of the right-most individual of the 
invading species, h y (t), for each row y. The average po- 
sition h(t) = (l/L y )J2 y hy(t) is then recorded for each 

time step, from which we extract the velocity [as h(t) 
approaches a constant slope for late times] . We ran each 
simulation until the front reached the end of the system. 
The longitudinal system size L x has no particular impact 
on the system's time evolution, although it constrained 
the maximal length of our simulations. 

One can also observe [Fig. [I], that as the front prop- 
agates, it also "roughens", i.e., the typical size of the 
fluctuations about the mean front position is increasing, 
before it reaches the steady-state for a given transverse 
system size L y . This kinetic roughening phenomenon 

111 in 

will be discussed in detail in Sec. IV. 

III. INVASION AS PROPAGATION INTO AN 
UNSTABLE STATE 

A. Mean-field equations and the asymptotic front 
velocity 

Taking the master equation corresponding to transi- 
tion rates in Eq. (Q, and neglecting correlations be- 
tween densities at different sites, yields dynamics of the 
ensemble-averaged local densities pi(x,t)=(ni(x,t)). We 
obtain 

Pi{x,t + 1) - pt{x,t) = [1 - pi(x,t) -p 2 (x,t)] 
x ~t M x '>*) - (2) 

x'€im(x) 
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FIG. 1: Snapshots of the moving and roughening front 
in the early-time regime (t—200 MCSS) and in the steady 
state (£=20,000 MCSS) for ai=0.50, a 2 =0.70, /i=0.20, and 
L y =200. White represents empty sites, while blue and red 
correspond to sites occupied by the residents and the invaders, 
respectively. 

where i = 1,2. Taking the naive continuum limit of the 
above equations, one obtains the (coarse-grained) equa- 
tions of motion 

dtPi(x,t) = ^[l-pi(x,i)-p 2 (x,t)]V 2 Pl (x,t) 
+ ai [1 - pi(x, t) - P2 (x, t)] pi(x,t) 
- PPt{*,t) , (3) 

i = 1,2. The spatially homogeneous solutions of these 
equations, (p*,/^), are (0,0), (1 — p/ai,0), and (0,1 — 
p/a.2). In the parameter regime of our interest, p < a% < 
«2j only the last solution (0, l—p/a.2) is stable. Thus, the 
propagation of a front separating the stable (0, 1 — pj 0.2) 
(invader dominated) and unstable (1 — p/a%,0) (resident 
dominated) regions amounts to propagation into an un- 
stable state Hf H, |2(| |2jJ |28j], a phenomenon that has 
generated a vast amount of literature [2£j since the origi- 
nal papers by Fisher Q and Kolmogorov et al. [B| . At the 
level of the mean-field equations, the front is "pulled" by 
the leading edge into the unstable state. Then, for a suffi- 
ciently sharp initial density profile [lll29||. the asymptotic 
velocity is determined by the infinitesimally small density 
of invaders that intrude into the linearly unstable region 
dominated by the resident species. Linearizing Eqs. © 
about the unstable state, pi = 1 — pja.\ + 4>i, P2 =02, 
one immediately obtains for the density of invaders 

d t M^t) = ^V 2 2 (x,t) + M f- 2 - - l) 2 (x,t) . (4) 
4 at \ai J 

Performing standard analysis 0, I27I l28l |29| on the 
above equation, we obtain the asymptotic velocity of the 



"marginally stable" invading fronts 

v* = —y / a 2 (a 2 - eei). (5) 

Oil 

The velocity v* above is the minimum velocity of a phys- 
ically allowed travelling wave, permitted by Eqs. Q, and 
is actually realized by deterministic nonlinear reaction- 
diffusion d yna mics for sufficiently sharp initial profiles 
I27], I28L l29l]. For further comparisons, we also note 
that the above asymptotic front velocity is approached 
as v(t) — v* — 0(1 /t) 29]. Also, as can be seen from 
Eq. JSJ , for small differences in the reproduction rates of 
the two species (a parameter of ecological significance), 
the front velocity scales as w*~(«2 — ot\) with = 1/2. 

It is important to note that the front velocity given by 
Eq. (JSJ, obtained by linearizing Eq. ©, fully reproduces 
the velocity obtained by numerically iterating the non- 
linear continuous-density mean-field equations [Eq. (|2*|)]. 
as can be seen in Fig. Pfa). This is a generic and pow- 
erful feature of deterministic pulled fronts where, despite 
the non-linearity of the full dynamics, velocities are com- 
pletely determined by the leading edge 0, |23 . 

B. Monte Carlo results for the front velocity 

We now present results for the discrete individual- 
based stochastic model defined by the transition rates 
Eq. (JIJ . We found that the front velocity is much smaller 
than that of the mean-field approximation as shown in 
Fig- E^b) [cf. Fig. |2t a )] ■ Further, for small differences in 
the reproduction rates, the velocity is described by an ef- 
fective power law w*~(a 2 - «i) e with 9 > 0.6 [Fig.^b)], 
slightly, but distinctly, greater than the exponent 6>=0.5 
of the mean- field case. Similar deviations from the mean- 
field exponent have been found in two-dimensional epi- 
demic [3(j and reaction-diffusion models [3l]]. The dis- 
creteness of the individuals [32l l33l l34| and noise 
have been shown to contribute to a velocity that is differ- 
ent from the mean-field description. The general belief 
is that fronts in stochastic individual-based (or particle) 
models are "pushed" , in the sense that the front veloc- 
ity is determined by the full non-linear front regi on, in- 
stead of an infinitesimally small leading edge [29||. a be- 
havior predicted by the mean-field approximation. Our 
two-species model provides an example for this generic 
behavior. 

An interesting feature of the invasion fronts is that 
their propagation velocity approaches the asymptotic 
value rather slowly (e.g. as 0(l/t) in mean-field mod- 
els [2^]). Thus, from application point of view, it is im- 
portant to establish how the front reaches its asymptotic 
velocity. In the next section, we are going to analyze both 
the temporal and finite-size approach toward the asymp- 
totic front velocity, along with other universal character- 
istics (such interface roughening |15l Hq^) of the model, 
using the framework of scaling in non-equilibrium inter- 
faces. 
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FIG. 2: (Color online) (a) Asymptotic mean-field front ve- 
locities for three different mortality rates for fixed ai =0.70, 
as a function of the difference of propagation rates, a.%— ot\. 
Symbols are obtained by numerically iterating the nonlinear 
mean-field equations of motion [Eq. J5J]. The solid curves 
running through the data points are the analytically obtained 
velocities in the "leading edge" approximation [Eq. 0], as 
described in the text. The straight dashed line indicates 
the slope 0=0.5, corresponding to the exponent of an effec- 
tive power law for small differences between the reproduction 
rates, (b) Asymptotic front velocities from Monte Carlo sim- 
ulations for the individual-based stochastic model for fixed 
a2=0.70 as in (a), with L y — 100, for three values of fi. The 
solid straight lines are the best-fit effective power laws in the 
region where the difference between the reproduction rates is 
small, corresponding to 6 = 0.61 ± 0.04, = 0.66 ± 0.04, and 
6 = 0.60 ± 0.04 for /i=0.1, /i=0.2, and /i=0.3, respectively. 
For comparison, the straight dashed line corresponds to the 
effective power law of the mean-field case with an exponent 
0=0.5. 



IV. INTERFACE ROUGHENING AND 
CORRECTIONS TO THE ASYMPTOTIC FRONT 
VELOCITY 

A. Dynamic scaling 

To extract the scaling properties of the roughening of 
the front, as can be seen in Fig. ^ we analyze the width 



of the front 



(w 2 (L yi t)) = 




where h y (t) is defined as the location of the leading indi- 
vidual in row y for the invading species. In what follows, 
we define L = L y , and investigate the scaling properties 
of the width (w 2 (L, t)} within the standard framework of 
dynamic scaling of non-equilibrium interfaces [TH Hf^ | . 

The temporal and system-size scaling of the width 
(w 2 (L,t)} typically identifies the universality class of 
the growing front. In finite systems, the width grows 
as {w 2 (L,t))^t 213 from early to intermediate times. At 
a system-size-dependent crossover time, t x ^L z , it sat- 
urates (reaches steady state) and scales as (w 2 at ) = 
(w 2 (L, oo)) ~ L 2a , where L is the transverse linear sys- 
tem size, a, /?, and z are referred to as the roughness, 
growth, and dynamic exponents, respectively. Further, 
these exponents are not independent, but obey the scal- 
ing relation a = [3z. The above temporal and system- 
size behavior, with the appropriate crossover time, can 
be captured by the Family- Vicsek scaling form |37| 

(w 2 (L,t)) = L 2a f(t/L z ). (7) 

For small values of its argument, f(x) behaves as a power 
law, while for large arguments it approaches a constant 

,/ n f x 2f} for i«l , . 

J[X> \ const, for x » 1 ' 1 ' 

yielding the scaling behavior of the width in the growth 
and steady-state regime, provided the scaling relation for 
the exponents, a = f3z, holds. 

As as shown in Figs. |3 our results show reasonable 
agreement with the exponents of the well-known KPZ 
universality class 03=1/3, and a=l/2) [HE El- Fa- 
ther, the scaled width, (w 2 (L , t)) / L 2a vs. t/L z produces 
good data collapse, as suggested by Eq. J7|), and confirms 
dynamic scaling for the invasion process [Fig.|3{a) inset]. 

B. Steady-state width distribution 

For further analysis, we also constructed the full dis- 
tributions P(w 2 ) of the steady-state width for different 
system sizes [i.e., the normalized histograms of the width 
obtained from the steady-state time series]. This ob- 
servable typically provides an additional strong signa- 
ture of the underlying universality class of the fluctu- 
ating and growing interface 38]. In particular, for the 
one-dimensional KPZ class, it has been obtained analyt- 
ically |38[ and can be written in the generic scaling form 
P{w 2 ) = {w 2 )- 1 ^>{w 2 /(w 2 )) with 

2 00 

$(a;) - \ VV-ir-Ve-^ 2 / 6 )" 2 *. (9) 
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FIG. 3: (Color online) (a) The average width as a function 
of time (on log-log scales) for various system sizes, averaged 
over 20 independent realizations, for Qi = 0.70, 0:2 = 0.80, 
and /x = 0.10. The dashed line corresponds to the one- 
dimensional KPZ power law with the exponent 2/3 = 2/3. 
The inset shows the scaled plot, (w 2 (L,t))/L 2a vs. t/L z us- 
ing the one-dimensional KPZ exponents, (b) Steady-state 
width as a function of the system size, L{= L y ), for the same 
clonal propagation and mortality rates as in (a). The solid 
line corresponds to the best-fit power law with the exponent 
2a = 0.95 ±0.01. 
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FIG. 4: (Color online) (a) Steady-state width distributions 
for ai=0.70, 0:2=0. 80, and /i=0.10 for various system sizes L. 
The solid curve is the scaled analytic width distribution of the 
one-dimensional KPZ class |3^|. Eq. ©. For comparison, the 
two-dimensional scaled KPZ width distribution is also shown, 
as obtained in Ref. |40j| (b) The same as (a) but on log-normal 
scales. 



stead of one transverse dimension in our model) . We give 
more discussion on this topic aspect in Sec. V. 



In Fig. 0] we show the scaled width distribution 
P(w 2 ){w 2 ) vs. w 2 1 ' (w 2 ) for various system sizes and com- 
pare it with the above analytic scaling function <&(x). 
Our data, again, strongly suggest that propagating pla- 
nar fronts in our two-species invasion model belong to 
the one-dimensional KPZ universality class. The large 
deviations and data scattered around peaks of the dis- 
tributions are due to sampling error. Steady-state time 
series for larger systems become strongly correlated, and 
so require excessively long simulations to generate suf- 
ficiently large statistically independent samples. Fig. 0] 
also shows the KPZ width distribution in two transverse 
dimensions 01 , offering a comparison suggested by a re- 
cent, somewhat counterintuitive conjecture |39| . That 
conjecture suggests that fronts which are pulled in the 
mean- field limit, belong to the one higher dimensional 
KPZ class than one would normally expect (i.e., two in- 



C. Temporal and finite-size corrections to the 
asymptotic front velocity 

Although the actual value of the velocity of the invad- 
ing front is not universal, Krug and Meakin showed 
that the forms of the temporal and finite-size corrections 
of the front velocity, v(L,t) = dh(L,t)/dt, are universal. 
More specifically, corrections to the asymptotic front ve- 
locity v* are given by jZJ 

(T ,s J v* - Cl t- 2 d-«)A for t^U- 
V ^ = \v*-c 2 L-^) for t>L* ' (10) 

where a and z are the roughness and dynamic exponents 
characterizing the universality class of the model, and C\ 
and C2 are non-universal constants depending on the pa- 
rameters and microscopic details of the specific model. 
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FIG. 5: (Color online) (a) Temporal corrections to the asymp- 
totic front velocity for L—S00 for different values of 02, with 
ai=0.50 and /i=0.20. The horizontal (time) axis is scaled as 
t~ 2 ^ 3 , motivated by the form of the corrections of the one- 
dimensional KPZ class [Eq. The inset enlarges the 
regime of (a) where the universal temporal corrections fol- 
low the KPZ behavior. The straight solid lines correspond 
to the linear scaling as a function of t~ 2 ^ 3 . (b) Finite-size 
corrections to the asymptotic velocity as a function of L~ , 
motivated by the universal corrections of the one-dimensional 
KPZ class, for ai=0.50, a 2 =0.70, and /i=0.20. The straight 
dashed line corresponds to the linear behavior as a function 
of 1/L. 



Translating their results to the easily measurable quan- 
tity h(t)/t (by integrating the above equation with re- 
spect to time and dividing the result by t), one also has 



h(L,t) 



t 



-2(l- a )/z for t < L z 

-c 2 L- 2(1 - a) for i>L z 



(11) 



where c[ = c\zj(z + 2a — 2). In particular, for the one- 
dimensional KPZ class, 2(1 — a)/ z = 2/3 and 2(1 — a) = 
1. Thus, the early-time temporal, and late-time finite- 
size corrections scale as (D(t~ 2 / 3 ) and C(i _1 ), respec- 
tively Our data for the velocity of the propagating front 
in the two species-invasion model follows these correc- 
tions very closely [Fig. [5], offering additional evidence 
that the front belongs to the one-dimensional KPZ class. 
Also, note that the temporal relaxation of the front ve- 



locity is in contrast to mean-field results where pulled 
fronts exhibit 0(t~ r ) corrections to the asymptotic ve- 
locity [23. 



V. DISCUSSION AND SUMMARY 

We studied front propagation in a two-species model 
for ecological invasion with preemptive competition, ap- 
plicable to clonal plants. We performed dynamic Monte 
Carlo simulations using the local transition rules [Eq. QJ] 
for initially flat linear fronts. We found the front veloc- 
ity significantly smaller than that of the mean-field ap- 
proximation. Also, for small differences between local 
reproduction rates, the asymptotic front velocity scales 

as an effective power law v*^(a2 — oli) with 9 > 0.6, 
an exponent slightly, but noticeably different from the 
mean-field value 9 = 0.5. The discreteness of the in- 
dividuals in our lattice model (or, the effec- 
tive density cutoffs in a continuum description) and noise 
[35l l36l| have been shown to contribute to the deviations 
from the mean-field results. More specifically, fronts in 
stochastic individual-based models, which in the mean- 
field limit converge to a pulled front solution, are instead 
"pushed" . Therefore, the front velocity is determined 
by the entire non-linear front region, instead of the in- 
finitesimally small leading edge alone [29| . Our model is 
an example for this type of behavior. 

We also investigated the universal features of the 
roughening of the propagating front, separating the in- 
vaders and the residents. We found that the front rough- 
ening in our two-species invasion model belongs to the 
one-dimensional KPZ class. We must also place our re- 
sults in the context of a recent conjecture by Tripathy et 
al. [39l | that propagating fronts, which in the mean-field 
limit are "pulled", exhibit KPZ scaling on a (d±+l) di- 
mensional "substrate" (where d± is the dimension of the 
space transverse to the direction of propagation), as op- 
posed to the naive expectation of a dj_-dimensional KPZ 
growth (i.e., two-dimensional instead of one-dimensional 
in our case.) That conjecture (3^ was based on field- 
theoretical arguments, but were subsequently shown by 
Moro 31] to apply only to systems where stochastic ef- 
fects are due to external fluctuations (such as, fluctua- 
tions in the parameters of the model). Most recently, 
it was argued [2^, and demonstrated [3l| that, in 
fact, fronts which in the mean-field limit are "pulled" , in 
the presence of internal fluctuations (i.e., systems with 
stochastic particle dynamics), belong to the "usual" d±- 
dimensional KPZ universality class. While corrections- 
to-scaling effects and system size-limitations can often 
hinder a high-precision determination of the exponents 
associated with interface roughening, the width distri- 
butions provide a very strong signature and aid in deter- 
mining the universality class of the interface [38| . To that 
end, we included the two-dimensional width distribution 
|40j| in Figs. supporting the conclusion that our two- 
species invasion model with stochastic particle dynam- 
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ics exhibits "standard" (one-dimensional) KPZ roughen- 
ing, in agreement with the most recent generic arguments 

As noted in the Introduction, although linear fronts 
are somewhat artificial in the context of ecological in- 
vasion, they do offer insight into more realistic scenar- 
ios. Consider, for example, that an advantageous allele 
or a competitively superior species is introduced through 
mutation within [jj, ITsj or through geographic disper- 
sal to j|| |23| a resident population, respectively. Intro- 
ductions occur rarely and stochastically in both space 
and time. Thus, small clusters of an advantageous allele 
or superior species can randomly "nucleate" and subse- 
quently grow. We have shown [aEJEJEll that the time 
evolution of such systems can be well described within 
the framework of homogeneous nucleation and growth. 
The growing clusters, on average, have radial symmetry 
and reach an asymptotic velocity v* for sufficiently long 
times. The corrections to the asymptotic radial veloc- 
ity of these circular fronts have two contributions: First, 
the typical length of the perimeter of the cluster scales as 
Lit) ~ 2nR(t) ~ t, where R(t) is the radius of the cluster. 



Since z > 1, t <C L z (t) for late times, i.e., the relevant 
KPZ correction for radially growing clusters is always of 
0(i -2 / 3 ) |4l|. Second, for long times, when the radius of 
the cluster is sufficiently large, the curvature introduces 
an additional ©(i?^ 1 ) ~ 0(t^ r ) correction, subdominant 
to the above KPZ correction. Thus, circular fronts are 
expected to reach the same asymptotic velocity as their 
linear counterpart, with the same leading-order 0(t~ 2 / 3 ) 
corrections for late times. 
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